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Abstract: We calculate the transition line of the first-order melting of vortex lattice in a 
three-dimensional type-II superconductor in fields of several Tesla, using the results from the 
density-functional theory of vortex melting in two- dimensions and a self-consistent Hartree 
treatment of correlations along the field. The result is in quantitative agreement with ex- 
periment. The temperature width of the hysteresis, the latent heat, the Debye- Waller factor 
and the magnetization at the transition are discussed. 
PACS numbers: 74.60. Ge 



Several recent experimental studies focused on the superconducting transition in un- 
twinned YBCO samples in fields of several tesla |jl|, |]. In contrast to earlier experiments 
in samples where the physics of the transition was dominated by inhomogeneities and where 
the transition appeared to be second order, the observed sharp hysteretic drop in resistivity 
in these very clean, strongly type-II materials in moderate magnetic fields suggests that the 
true vortex melting transition might be first order. This is somewhat surprising when one 
recalls that the mean-field Abrikosov theory predicts a second order phase transition for a 
homogeneous type-II superconductor in magnetic field. Thus here one encounters another 
possible example of strong thermal fluctuations changing the order of transition. An early 
suggestion that this might happen in a type-II superconductor in vicinity of Hc2{T) came 
from the renormalization group analysis close to the upper critical dimension d^p = 6, as 
well as from the analysis of the theory with infinite number of order parameter components 
in 4 < d < 6 [0. This is inadequate however for the physical three-dimensional (3D) samples 
which are below the lower critical dimension diow = 4 in the problem considered there. Usu- 
ally, the transition in the vortex system is described by the harmonic theory of vortex lattice 
and by invoking phenomenological Lindemman criterion to locate the melting point [0, 
This is however a suspect starting point if one is interested in describing strong fluctuations 
near Hc2(T) and is more appropriate for the low-field or T region of the phase diagram. 
As emphasized by Moore ||^, in this description one starts from the Abrikosov lattice solu- 
tion for the order parameter which is unstable with respect to harmonic shear modes of the 
lattice at any finite temperature in two and three dimensions. It is therefore inconsistent to 
simply assume this ordered low-temperature state. Furthermore, the numerical constant in 
Lindemman criterion needs to be chosen phenomenologically, and the requisite number can 
actually differ by orders of magnitude from one material to another \W\. 

Recently, a novel approach to the problem has been formulated that encompasses the 
difficulties mentioned above by relying on a new physical picture of the phase transition in 
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strongly type-II system in high field |Tl|. Unlike the renormalization group study of ref. 5 
where the superconducting Abrikosov transition arises from the growing phase correlations 
in the directions along the magnetic field, here the positional correlations orthogonal to the 
field drive the transition. The right paradigm is the two-dimensional (2D) problem. In high 
magnetic field (to be specified later), the original Ginzburg-Landau partition function in 
the symmetric gauge is equivalent to a system of classical particles interacting with long- 
range, multi-body forces |]1T|. It is a scale invariant, incompressible system which undergoes 



a weak first-order freezing transition, not unlike one-component Coulomb plasma in 2D. 
The point of transition, the latent heat and other relevant quantities at the transition can 
be quite accurately determined by using the density-functional theory of solidification [|12] . 



The low-temperature phase is not the familiar Abrikosov vortex lattice but a charge-density- 
wave (CDW) of Cooper pairs, with a weak periodic modulation of Cooper pair density but 
with no long-range phase coherence [0, |13|. In this paper we extend the theory to 3D 
superconductor and use it to calculate the phase boundary in H — T phase diagram. We find 
a very good quantitative agreement with experiment. Also, we make a number of predictions 
about the typical physical quantities at the transition and analyze the temperature width of 
the hysteresis as a function of the field. 

Consider the Ginzburg-Landau partition function for an anisotropic, homogeneous, 3D 
superconductor in strong magnetic field, with fiuctuations of the magnetic field neglected 
(Ginzburg-Landau parameter k » 1). For the fields above Hb ~ {9/16)Hc2{0)iT/Tco), 
where 9 = A^jf/Ax (Ae// = 2Xab{0)'^/d, At = 0g/167r^Tco) is the Ginzburg fiuctuation 
parameter the fiuctuation spectrum is dominated by the lowest Landau level modes (LLL) 



14|. In YBCO 6 ~ 0.04 and Hc2{0) ~ 150T so the approximation which retains only the 
LLL modes should be appropriate for the fields larger than ~ 0.3T. The partition function 

is Z = J D[ip* ,ip]exp{—S), and 

S=^ld'rdz[^\d,i;\' + a'{T,HM' + ^\^P\'] (1) 



where a'(T) = a(T){l — H/Hc2{T)), and a(T) = a(T — Tco), /5 and 7 are phenomenological 
parameters. It is convenient to rescale the fields and lengths as {2dp2Trl'^ /TY^^ip ip, 
r/(l\^27i) — i> r and z/d z, where / is the magnetic length for charge 2e, and c? is a 
microscopic length along the field (typically the spacing between the pairs of CuO planes). 
The action in the exponent then becomes 

S = J d'fdz[g,\d.^\' + gM' + \\i^\% (2) 

with ip which is restricted to the LLL, and the whole thermodynamics is determined by two 
dimensionless coupling constants g^/a = {l/d"^, a} x ^rrPd/ PT. In the 2D = 0, and 

the system described by the partition function (2) undergoes a weak first-order transition 
into CDW phase of Cooper pairs at ga = gM = —6.5, as seen in Monte-Carlo simulations 
n|, and found in density-functional theory |T2|. The CDW has a triangular 



modulation with the period set by the magnetic length: a = lyAn/y/S. Small deviations 
from this periodicity are expected in principle, but do not change our results substantially 
and will not be considered here. We assume that the vortex transition in 3D system is 
driven by the same mechanism of growing positional correlations between vortices. In 2D, 
ga oc (i^/^, where d is the film thickness. In 3D, this length is replaced by some temperature 
and field-dependent length A over which the vortices are "straight" along the field direction. 
A provides a short length scale along the 2;-axis, just like the magnetic length does in the 
X — y plane. Here we assume that this length is not very different from the superconducting 
correlation length along the field, ^||. Thus we make an ansatz that there is a single coupling 
constant which describes the physics of the transition in the 3D regime (i. e., for ,^[[ > 1 in 
units of d): gaC,^/'^- The phase boundary is then determined by the condition 

9m = 5'aC||''^ (3) 

for > 1. To explicitly determine the transition line in the H — T phase diagram we use the 
correlation length as obtained in the self-consistent Hartree treatment of the theory defined 



by eq. 2: 



A, = ( sl/2 U) 

where the thermal average appearing is determined by the equation 

^ 1"^' 2^ /-oo + (^,+ < |^|2 > /4- 
The expressions 4 and 5 look the same as they would in a purely ID theory. This is a 
consequence of the strong magnetic field, which causes D ^ D — 2 dimensional reduction 
in the theory on the Hartree level. The full theory as given in eq. 2 is, of course, not 
ID, due to the non-local constraint on the fields to be entirely in the LLL. In fact, the 
transition we are describing comes exactly form those lateral, intra-LLL correlations which 
would be undetectable in a simple Hartree theory. ^||, however, is expected to be reasonably 
well described by the Hartree theory, since the dimensional reduction is exact for 
correlator [0, [T^ . 

From the melting condition 3 and using the expressions 4 and 5 one obtains the equation 
for the transition line 

t + + (^^^y^)^/ = 1, (6) 
where we rescaled the field and the temperature as t = T/Tco, h = H/Hc2{0), and the 



constant c is determined by the numerical value of qm as = glfiyl + l/f?!/ ~ = 10.5. 
The equation for the transition line 6 is our main result. It depends only on a particular 
combination of the material constants and on a pure number qm- Its accuracy is restricted 
to the region where ^|| > 1 since otherwise the layers would decouple and the description in 
terms of 2D melting would be appropriate. In the region of 3D to 2D crossover better results 
would be obtained by using the Lawrence-Doniach model. We chose anisotropic Ginzburg- 
Landau partition function instead because the transition line can be found in simple closed 
form, which still is accurate over a large portion of the phase diagram. 

In figure 1. we compare the calculated transition line with the experimental results on 
superconducting transition in YBCO in the magnetic field parallel to the c-axis [0]. We set 



Tco = 90K, Hc2{0) = 148. 14T (as found by linear extrapolation from small fields), k = 50, 
^c(O) = 3A, and the best fit is obtained with ^a6(0)/^c(0) = 5.8 which agrees well with what is 



known about this material [T^ . The results of experiments in refs. 1 and 2 are essentially the 
same and are not shown. Note that the distance between CuO planes d = 12A completely 
cancels out in the determination of the transition line. According to our calculation, YBCO 
is three dimensional everywhere along the transition line in the range of the fields studied. 
One may also expect that the numerical value of qm should be renormalized from its 2D 
value that we used due to the difference between the lengths A and ,^||. This question can be 
settled by knowing more precisely the values of parameters k and anisotropy in eq. 6. The 
present calculation suggests however that this difference is not substantial. 

The density-functional theory predicts that the CDW phase can be superheated up to 
QsH = —6.25 in the 2D vortex system. If we neglect the supercooling of the vortex liquid 
upon lowering the temperature and take only the superheating of the solid phase into account 
we may calculate the thermodynamic hysteresis width in temperature as a function of the 
magnetic field for 3D transition. The result is shown on Figure 2. The functional dependence 
of the calculated hysteresis width agrees with the observation in resistivity measurements in 
refs. 2,3,4 but the result is roughly an order of magnitude larger. We attribute this to the 
effect of disorder in the experiment, which tends to reduce the width of hysteresis. However, 
for a given sample and a small current, the hysteresis widths in temperature at different fields, 
when multiplied by a suitably chosen fudge factor, agree well with our curve (see Figure 2.). 
The reason for this is the following: in a homogeneous 2D sample, the hysteresis width 
^thyst{h) oc Ag = gsH — Qm-, for small Ag. A weak point disorder cuts off the crystalline 
order at large but finite Larkin-Ovchinikov length ^lq- For a crystallite of that size both gM 
and gsH increase proportionally to ^^q, but their difference Ag decreases. Thus the primary 
effect of weak disorder is to decrease the numerical value of Ag from its thermodynamic value 



of 0.25 ||T2[, and therefore decrease Athyst{h). The functional dependence of the temperature 
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width of the hysteresis on the field however is not affected and hence the above behavior. 

Ahhough the transition under consideration is into a 3D phase of periodically modulated 
density of Cooper pairs, we expect the relevant quantities at the transition not to be very 
different from 2D case. This comes as a consequence of the same mechanism of the transition 
which is effectively 2D in nature. Thus, the latent heat per vortex and per layer should be 
~ O.S/cbTm [|12|. Debay- Waller factor ly^G) = |p(G')p at the transition should be given by 
p{G) = 0.72 exp(— A^G^) with A = 0.47, where G is a reciprocal lattice vector in units where 
1 = 1. The magnetization is expected to have a discontinuity of roughly 1% of its value at 
the transition [1^. The solid phase below the transition is expected to have at most a power 
law superconducting order |P, |ll| , but it will behave as a superconductor for most practical 
purposes-For example, its ohmic resistivity will be extremely low. This is the consequence 
of YBCO being effectively an anisotropic 3D superconductor in the region of fields and 
temperatures where the experimental data are currently available. This leads to an enhanced 
density modulation of the SCDW below the transition with the resulting large increase in 
the range of superconducting correlations. At even higher fields, where the layered structure 
of YBCO becomes more pronounced, we expect the range of superconducting correlations 
to become substantially shorter, perhaps leading to a clearly identifiable SCDW state with a 
long-range 3D positional order but only a finite resistivity. In this respect, a large variety of 
the HTS materials which are more anisotropic than YBCO appear particularly promising. 

Finally, we should stress that our results leave open a possibility that qmI^D) is sig- 
nificantly different from qmC^^D), with the corresponding modifications in the latent heat, 
characteristic vortex bending length, magnetization discontinuity, etc. The problem is that, 
while our theory successfully captures lateral correlations of the 'microscopic' GL-LLL the- 
ory which are responsible for the SCDW transition both in 2D and 3D, it still requires as an 
outside input the vortex structure factor near the transition. In 2D we were able to obtain 
this needed information from the numerical Monte Carlo simulations of the GL-LLL the- 
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ory. Unfortunately, such numerical results are not available at present for the 3D or layered 
GL-LLL model. 

In conclusion, we calculated the transition line of first order vortex lattice melting in 
strongly type-II superconductor in magnetic field of several Tesla. The agreement with 
the experimental results on untwinned YBCO is excellent in the whole range of fields of 
0.1 — lOT. The temperature width of the hysteresis, the latent heat, the Debye-Waller 
factor at the transition and the discontinuity of the magnetization are discussed. Our results 
strongly support the conclusion that the vortex-liquid to solid transition in homogeneous 
strongly type-II superconductor is first order. 

This work has been supported in part by the NSF Grant No. DMR-9415549. 
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Figure Captions: 

Figure 1. H-T phase diagram for YBCO (t = T/Tco). The full line represents the 
melting line calculated from eq. 6 with parameters described in the text. Crosses are the 
experimental points from ref. 3. The dashed line is the mean-field Hc'ziT) line. 

Figure 2. The field dependence of temperature width of the hysteresis due to superheat- 
ing of the low-temperature phase. The full line represents the calculated hysteresis width. 
Triangles and crosses are experimental points of ref. 2 for two samples with different degrees 
of disorder. The temperature widths are multiphed by factors 5 (triangles) and 10 (crosses) 
to fit the curve. 
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